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A dilute bose gas in a quasi two-dimensional harmonic trap and interact- 
ing with a repulsive two-body zero-range potential of fixed coupling constant 
is considered. Using the Thomas- Fermi method, it is shown to remain in the 
same uncondensed phase as the temperature is lowered. Its density profile 
and energy are identical to that of an ideal gas obeying the fractional exclu- 
sion statistics of Haldane. 
PACS: 03.75.Fi, 05.30.Jp, 67.40.Db, 05.30.-d 



I. INTRODUCTION 

The quantum and thermodynamic properties of atom traps with reduced dimensionality 
have been studied by several groups In a recent review Dalfovo et al. ^ note that 

the phase transition of bosons in quasi two-dimensional atom traps is an important issue 
for future investigations. In this paper, we consider bosons which are trapped in a two- 
dimensional harmonic oscillator and interact by a repulsive zero-range pair potential of a 
fixed strength. Using the finite-temperature Thomas- Fermi (TF) method, we find that there 
is no strict phase-transition for such a system, no matter how weak the repulsion. As the 
temperature is lowered, the chemical potential rises smoothly, and only at zero temperature 
does it match the lowest energy level of the trap. By analyzing the equation of state of the 
interacting bose system in the complex fugacity plane, we show that the system remains 
in a single phase at all (real) temperatures. We demonstrate that the area integral of the 
interaction potential in two dimensions is related to the statistical parameter of Haldane's 
fractional exclusion statistics (FES) ||^. In particular, the thermodynamic properties of the 
bose system interacting with a zero-range repulsive potential are found to be of the same 
form as that of an ideal FES gas. 

Consider bosons in an oblate three-dimensional trap with the oscillator frequencies 
= Uy = uj Uz. For a dilute Bose gas at low temperatures we take the usual three- 
dimensional delta function pseudo-potential with strength {Anh'^ a / m) , where a is the 
s-wave scattering length, and which in first-order Born-approximation yields the correct 
low-energy quantum scattering result. For low enough temperatures, and assuming that 
the expectation value of the pseudo-potential is much smaller than huJz, we may restrict 
the Hilbert space in the z-direction by setting the oscillator quantum number = 0. The 
quantum dynamics is then determined by an effectively two-dimensional Hamiltonian, given 
by 
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where the momenta and coordinates are planar vectors. Such a Hamihonian has been studied 
recently by several groups 0J^|1^ and also by Mullin ||TI| , p!2| . The dimensionless coupling 



constant g in the two-dimensional delta function interaction may be related to the scattering 
length a for s-wave scattering in three dimensions. This is done by taking the expectation 
value of the three-dimensional delta function with respect to the harmonic oscillator one- 
dimensional ground state wave function with = 0. We then obtain the effective two- 
dimensional interaction given in Eq. (j^), with the dimensionless coupling constant 




9 = \-^ ■ (2) 



Here bz = \Jh/ mojz is the length scale of confinement in the z-direction. This quasi two- 
dimensional description of the actual three-dimensional trap is valid only for a dilute system 



in which the scattering length a <^ bz <^ b where b = \Jh/muj is the oscillator length 
scale in the x-y plane. It follows from Eq. (0), therefore, that g <^ 1 m Eq. (0). For an 
estimate of the parameter g, take 6 = 3 micron, which corresponds to huj ~ 0.6 nK. An 
ideal gas in a two-dimensional harmonic trap would have BEC at the critical temperature 
T° = {Q/Ti'^f'^N^I'^huj . Taking = 10^ Rb^^ atoms, this corresponds to T° ^ 48 nK. 
We must have hu^ ^ ksT^ for the device to be effectively two-dimensional. Assuming 
ujz — 10^ a;, and the scattering length a for Rb^^ atoms to be 5.8 nm |]13[, the estimate of the 
interaction strength g from Eq. (§) is about 0.05. Thus the model Hamiltonian (|l]) is real- 
izable experimentally for small g, although we shall also present the results for = 1 to see 
how the quantum statistics is altered. We emphasize that the two-dimensional delta func- 
tion interaction in Eq. (|l|) is a dimensionally reduced form of the three-dimensional pseudo- 
potential, and is not to be regarded as a pseudo-potential derived from two-dimensional 
scattering theory. Later, we shall need to take the thermodynamic limit of this quasi two- 
dimensional system to investigate the question of strict phase transition. This is done by 
keeping uOz at a fixed appropriate value and letting A^ oo, — 0, keeping N^uo constant. 
Note, however, that in the alternative scheme of keeping the ratio uJz/^ fixed at a large 
value, and then taking the thermodynamic limit, would result in a three-dimensional geom- 
etry. This is because the condition fiujz 3> ksT^ would not be satisfied. We now proceed to 
do a Thomas-Fermi calculation for this two-dimensional system. 



II. THOMAS-FERMI CALCULATION 



In BEC, it is customary to write the density of bosons to be n(r) = no(r) +nj'(r), where 
no is the condensate density, and the density of particles occupying states other than 
the ground state 0. Consider temperatures above the critical T^. of the interacting system. 
Then there is no condensate, and n(r) = nT(r). In this situation, the one-body potential 



generated by the above zero-range interaction is ?7(n(r)) 



gn{r). Note that for T > T^, 



the more sophisticated Popov approximation also reduces to the TF approximation: 



n(r] 



exp[(£ + im^V2 + '-^gn{r) - - 1 



(3) 
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After performing the ]9-integration analytically, the number- density of particles {(3 = 
l/{kBT)) is given by 

^g-/3m..V2/2 _ 2e""sinh- , (4) 

where we have put u{t) = 27rh'^n{r)P/m, the fugacity Z = exp(/3/i), and a = {g — 1/2). The 
number of atoms in the trap is given by 

N = [ uMd^r . (5) 

To obtain n(r) for a fixed A^, Eqs. and have to be solved self-consistently. For g = 0, 
these cannot be satisfied for T < Tj^, and the occupancy of the lowest quantum state (in this 
approximation at zero energy) has to be taken into account explicitly. This means that the 
condensate density uq occupying the lowest level is macroscopic below T < T°. In Fig. 1, we 
show how this results in a discontinuity in fi at for g = 0. (In numerical work throughout, 
we use arbitrary units [a.u.] with h = m = 1 and set N = uj = 1.) 
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FIG. 1. Chemical potential /i versus temperature T /T^. Self-consistent solutions of the 
temperature-dependent TF equations (|^,§) are obtained for various values of the interaction 
strength g 



For a nonzero positive Eqs. @) and (^) can be self-consistently solved right down 
to T = (see Fig. 1). This contradicts MuUin's claim |T2[ that there is no self-consistent 



solution of the TF equations below a nonzero T^. Note that the T — * limit of Eq. (1) yields 
a nonzero spatial density only within the classical turning point tq: 

no{r) = ^J^2^ (/^ - ^mu^r^^ , r < Tq , (6) 

where tq = \f^^- This gives = fj? /2g{fiujY , so that the chemical potential at zero 
temperature is given by /i(T = 0) = y/2g N^^'^fiuj . Thus, the self-consistent solution has 



3 



led to the T TF result (Eq. (|D) for uq that one would obtain starting from the Gross- 
Pitaevskii density functional In passing we note that the zero-temperature density for 
a spin-less free Fermi gas in a two-dimensional harmonic potential is also given by Eq. 
with the proviso that the occupancy of each quantum state is 1/g instead of unity. Indeed, 
even though it is not experimentally feasible to extrapolate g to unity, it is amusing to note 
that putting g = 1 in Eq. @ would yield exactly the finite temperature noninteracting 
fermionic result 

TTl / X \ 

nir) = In [I + exp(/i - -mu;^r^)pj . (7) 

These results are not accidental, and their significance (specially for small g) will presently 
be discussed. It is also worth noting that the plots in our Fig. 1 are very similar to those of 
Fig. 2 for a noninteracting Haldane gas in Ref. [1^ . 



It is important to note that the situation is very different for an interacting gas in a 
three-dimensional trap, where the TF equation analogous to Eq. (0) has no self-consistent 
solution for T < Tc. One then has to solve the Bogoliubov equations self-consistently, 
say in the Popov approximation 0. For T > this procedure reduces to solving the 
temperature-dependent TF equation. Qualitatively, the calculated curves of yU vs T for the 
three-dimensional trap (shown in Fig. 26 of 0) look similar to our Fig. 1. This is somewhat 
deceptive, since the curves in our Fig. 1 are obtained in the TF approximation for all T, 
whereas the Popov approximation is used for T < Tc in the three-dimensional case. The 
behavior of the so-called "release energy" , which is the total minus the harmonic oscillator 
energy, has a very interesting structure below and at in the three-dimensional case. 
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FIG. 2. Release energy versus temperature T /T^ from the self-consistent solution of the 
temperature-dependent TF equation. For comparison, the noninteracting g = Q case is also 
shown 
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for Rb atoms, who find a dis- 
as is shown in Fig. 2 for 



This is also confirmed experimentally by Ensher et al. 
continuous behavior of the deduced specific heat. By contrast, 
the two-dimensional problem, no such structure is found for nonzero g in the calculated 
release energy. Experimentally, in three-dimensional traps, the appearance of a condensate 
is manifested as a sharp peak in the spatial density profile of the gas [|T^. For the interact- 
ing two-dimensional gas, as shown in Fig. 3, the central density rises gradually rather than 
abruptly with the lowering of temperature. (This, however, may be a limitation of the TF 
model, see 01.) 
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FIG. 3. Self-consistent TF density distribution, with unit normalization, as a function 
of the radial distance r for temperatures above and below T°. The noninteracting g = 
case is shown for comparison 



We now show analytically, starting from Eqs. (§J^), that for g > there can be no 
phase transition. Formally, these equations bear a striking similarity with Eqs. (29,30) of 
Sutherland's classic paper on the thermodynamics of the Calogero- Sutherland model 
(CSM) |T9|. We exploit this for our proof, going to the thermodynamic limit —>■ oo, 
a; ^ 0, with N^^'^uj = 1. In this limit, the r-dependence in Eq. (H) drops out, and we get 



Z = [exp{gu) - exp{{g - l)u)] 



(8) 



From this, we see that for g = 0, u = — ln(l — Z), so u has a branch-point at Z = 1, i.e., 
at /i = 0. This is where BEC takes place for the ideal Bose gas. For g = 1, on the other 
hand, Eq. (^ gives u = ln(l -|- Z), with a branch point at Z = —1. Similarly, for g = 1/2, it 
is easily checked that u has branch points at Z = ±2i. More generally, the location of the 
branch points of u (which is proportional to the spatial density) in the complex Z-plane as 
a function of g are shown in Fig. 4, which is identical to Sutherland's Fig. 1 for CSM in [jl8| . 
The only singularity on the positive real axis is at = 0, and for g > there is no phase 
transition. 
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FIG. 4. Location of the branch points of u in the complex Z-plane as a function of g 
id. @) 



III. EQUIVALENCE TO FRACTIONAL EXCLUSION STATISTICS 

The two-dimensional delta function interaction considered in Eq. (|l|) is very special in the 
sense that it is scale-independent A two-body potential that varies as the inverse square 
of the distance between the two particles also shares this property. In one dimension, the 
latter is the much-studied exactly solvable A^-particle CSM, which is known to obey [^0| -|22 



Haldane's fractional exclusion statistics (FES). By this we mean that either one can solve 
the A^-body interacting bosonic or fermionic problem with the pair-wise inverse square in- 
teraction, or map it onto a noninteracting set of A^-particles obeying FES. The statistical 
occupancy factor for an ideal FES gas is given by p2|,^ rji = {wi + g)^^, where 



wf{l+w,f-^ = exp[{e,- ^J)P]. (9) 

We note that FES is characterized by a "statistical factor" g, which will presently be identi- 
fied with the dimensionless interaction strength of the delta function interaction in Eq. (P. 
The occupancy factor in FES reduces to the usual bosonic or the fermionic one for (7 = or 
1, respectively. 

In order to demonstrate explicitly the equivalence between our two-dimensional interact- 
ing boson model and FES, let us consider the analogue of Eq. (|]) for the number density of 
a gas of non-interacting "haldons" , 

where w depends on p via Eq. . The number density can be rewritten as an integral over 
w by noting that Eq. (|) gives \glw + (1 — + vS)\dw = {P/2TTm)d'^p, so that 

m f°° dw m 1 + Wo 
n r = T— / — ; = ^ m , (11) 

^ ' 27ih^p Jwo w{l + w) 27ih^P Wo ^ ^ 
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with wq = w{p = 0) given by 



(12) 



Eliminating Wq from Eqs. ([Tl|) and (^) just gives Eq. @, which shows that the number 
densities of the two models are the same. We proceed to prove that the same is the case for 
the energy density, which for the FES gas is given by 



p mtu 

— H 

2m 2 



w + g 

Changing variables from p to w as above and using Eq. (|^) to substitute 

mco'^r^/2] = Pn + glnw + {1 — g) ln(l + w), 
the energy density takes the form 

dw 



E{t) 



m 



gin 



w 



1 + w 



+ ln(l + w)+pfi 



2T^h^P^ Jwo W{l + W) 

The first term of the integrand is a total derivative, giving a contribution 

2 " dw 



(13) 



(14) 



(15) 



-m^/47rr/?^) In 



Wo 



1+Wo 



-nh g/m){n{r))^ 
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2/32 



wo 



W[l + W 



-n(r) (16) 



to the energy density. With the variable substitution w = Wq + {1 + Wo){exp{y) — 1), the 
remainder of the integral (15) takes the form 

nfi^g , , ..2 f°° dy 



{n{r)y 



m ■ " 2nh^p^Jo {l + Wo)ey-l 
Here, Wq can be eliminated by using the relation 

1 + Wo = exp |/3 muj'^r^/2 + {2TTfi^ g /m)n{r) — | 



[y + ln(l + Wo) + PfA . 



(17) 



(18) 



which is found by combining Eqs. (|Tl|) and (^. Further identifying y = (3p'^/2m and 
collecting all terms, we find the final form of the energy density integral. 



Eiv) 



d^p/{27rny 


i£ + im^V2 + 






exp[(|i + |ma;V2 + '-^gnir] 


-/i)/3]- 


- 1" 



(19) 



But this is just the energy density of the interacting boson gas. Note that Eq. (|TBp brings in 
a term ~ gn{T) on the right-hand side of Eq. (0), which combines with the one in Eq. ( p!6[) 
to give the correct coefficient nfi^g/m in Eq. (|TP|). Thus, we have demonstrated that the 
interacting bosons may be mapped on to noninteracting particles obeying FES. Although 
above we have assumed a harmonic trap, note that our proof holds for any shape of the 
confining trap potential. 

For a dilute gas in the thermodynamic limit obeying FES in a D— dimensional system, 
the statistical factor g may be related ||2^;25] to the high-temperature limit of the second 
virial coefficient B2 [^ : 
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9-\ = 2^/'i?2 . (20) 

For the parameter g from Eq. (pO]) to be meaningful, it should be temperature- independent. 
This is possible if the interaction potential is scale-independent. This requirement is satis- 
fied by an inverse-square potential in any dimension. In one dimension, the corresponding 
potential is the CSM mentioned before. Interestingly, the two-dimensional case is very spe- 
cial, since in this case the area-integral of any potential, when it exists, turns out to be 
proportional to the statistical parameter. To see this, consider, for (D = 2), the /5 — > limit 
of B2 given in 



Tfi 1 

S, = ^M„±j, (21) 

where Mq = / dPrV{Y). We have assumed that the potential V is well-behaved, so that 
the expansion exp(— /5y) = [1 — (3V + {(3VY /2 — ...] is valid. Terms containing / (PrV"' 
for n > 2 do not contribute in the high-temperature limit, since Z?'-"^^^ 0. Considering 
bosons in D = 2, Eq. (^) gives Mo = {27ih'^/m)g . The zero-range effective potential in 
Eq. ([^) has precisely this moment Mq, and its strength g is thus the statistical parameter 
for Haldane statistics. 



IV. DISCUSSION 

In summary, we have shown, using the Thomas-Fermi method, that the effect of a repul- 
sive zero-range interaction between the atoms in the "two-dimensional" trap is more drastic 
than the corresponding three-dimensional case. Even in the thermodynamic limit, there 
is no strict phase transition with a repulsive zero-range two-body interaction, as shown by 
the behavior of the spatial density in the complex Z-plane. We have also demonstrated 
that such interacting bosons in a two dimensional trap have identical bulk thermodynamic 
properties as those of noninteracting "haldons" in the same trap, obeying FES. 

The above conclusions have been reached using a two-dimensional zero-range interaction 
with a fixed coupling constant, which was obtained by a dimensional reduction of the three- 
dimensional pseudo-potential, as explained in the introduction. It should be noted, however, 
that in a strict two-dimensional problem, the pseudo-potential is more complicated. As 
emphasized by Schick , the scattering cross section for binary collision of two bosons due 
to a hard-disc potential of radius a diverges in the low-energy limit : 



(ka) « 1 , (22) 



where k is the relative wave number. This is a very different behaviour from the three- 
dimensional problem, where the total cross section goes to a constant value proportional 
to a^. Thus, in strictly two- dimensions, the strength of the pseudo-potential should be 
/c-dependent. 

Another way to define a pseudo potential for a very dilute gas would be to have a 
zero-range form with a density-dependent strength such that in Born approximation it re- 
produces the ground-state energy per particle to the lowest order. Schick had shown that 
this is given by E/N = — 27m(ln na^)~^, where n is the average number density of the bosons. 
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Shevchenko in a detailed analysis of bosons in a trap, thus adopted a zero-range po- 
tential with the g of our Eq.(§) replaced by the density- dependent factor ln~^(l/(na^)). He 
came to the conclusion that whereas this repulsive interaction prevents strict BEC, super- 
fluidity sets in at a temperature close to the critical value T° = (6/7r^)^/^A^^/^^a; of 
the ideal gas. This transition to superfluidity for the compressed fraction of the gas near 
the center of the trap takes place through the Berezinskii-Kosterlitz-Thouless (BKT) |]30[| 
phenomenon of bound vortex pairs, much as in the translationally invariant two-dimensional 
interacting gas. 

In our simple Thomas-Fermi treatment of the two-dimensional problem with a zero- 
range pseudo potential of constant coupling strength, we cannot examine the fluctuations in 
the phase of the "quasi-condensate" nor can we comment on the formation of BKT- 



vortices. Despite these limitations, we find, for our model, the transformation of Bose to 
Haldane statistics through the zero-range interaction. This, arguably, is the most interesting 
point of our paper. 
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